clc;
clear;
B=[2,6,-1,-2;
    5,-1,2,3;
    -3,-4,1,10;
    5,-2,-3,8];

A=[-4,-6,-4,-1;
    1,0,0,0;
    0,1,0,0;
    0,0,1,0];
[V,D]=eig(A,B),
norm(A*V-B*V*D)